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Heavy ion reactions and other collective dynamical processes are frequently described by different 
theoretical approaches for the different stages of the process, like initial equilibration stage, inter- 
mediate locally equilibrated fluid dynamical stage and final freeze-out stage. For the last stage the 
best known is the Cooper-Frye description used to generate the phase space distribution of emitted, 
non- interacting, particles from a fluid dynamical expansion/explosion, assuming a final ideal gas 
distribution, or (less frequently) an out of equilibrium distribution. In this work we do not want to 
replace the Cooper-Frye description, rather clarify the ways how to use it and how to choose the 
parameters of the distribution, eventually how to choose the form of the phase space distribution 
used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with 
the freeze-out problem, while the discussion of transition between different stages of the collision is 
applicable to other transitions also. 

More recently hadronization and molecular dynamics models are matched to the end of a fluid 
dynamical stage to describe hadronization and freeze-out. The stages of the model description can 
be matched to each other on spacetime hypersurfaces (just like through the frequently used freeze- 
out hypersurface) . This work presents a generalized description of how to match the stages of the 
description of a reaction to each other, extending the methodology used at freeze-out, in simple 
covariant form which is easily applicable in its simplest version for most applications. 



I. INTRODUCTION 

Relativistic heavy ion reactions exhibit dominant col- 
lective flow behaviour, especially at higher energies where 
the number of involved particles, including quarks and 
gluons, increases dramatically. At intermediate stages 
approximate local equilibrium is reached, while the ini- 
tial and final stages may be far out of local equilibrium. 
Also, different stages may have different forms or phases 
of matter, especially when Quark Gluon Plasma (QGP) 
is formed. 

The need to describe and match different stages of a re- 
action was realized by the development of the final freeze- 
out (FO) description in Landau's fluid dynamical (FD) 
model [lj. Then it was improved by Milekhin and 
a covariant simple model was given by Cooper and Frye 
In all these models the FO happened when the fluid 
crossed a hypersurface in the spacetime. 

At early relativistic heavy ion collisions, the initial 
compression and thermal excitation was described by a 
compression shock in nuclear matter. This was already 
pointed out by the first publications of W. Greiner and 
E. Teller and their colleagues d, and the shock took 
place crossing a spacetime hypersurface (e.g. a relatively 
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thin layer resulting in a Mach cone) . When sudden large 
changes happen across a spacetime front the conservation 
laws and the requirement of increasing entropy should be 
satisfied: 

[N'dvp] = ; (1) 
[T^dcTfj] = ; (2) 
[S^da^] > (3) 

where — nu* 1 is the baryon current, = su M is the 
entropy current, T 1 *" is the energy momentum tensor, 
which, for a perfect fluid, is given by 

= (e + P)u»u v - Pg**" , (4) 

where e is the energy density, P is the pressure, s is the 
entropy density, and n is the baryon density of matter. 
These are invariant scalars. The da^ is the normal vector 
of the transition hypersurface, is the particle four ve- 
locity u^ 1 — 7 (1, v x , v y , v z ) = 7 (1, v), normalized to +1. 
The square bracket means [a] = a\ — an, the difference of 
quantity a over the two sides of the hypersurface. The 
metric tensor is defined as g^ v = (1,-1,— 1,-1). We will 
also use the following notations: w = e + P, j = iV^doy, 
is the invariant scalar baryon current across the front, 
X = (e + P)/n 2 = w/n 2 is the generalized specific vol- 
ume, v 2 = v 2 — v 2 x + v 2 + v 2 z , and fx, v... = 0, 1, 2, 3, 
i,3— = 1,2,3 = x,y,z. 
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For a perfect fluid local equilibrium is assumed, thus 
the fluid can be characterized by an Equation of State 
(EoS), P = P(e,n). Eqs. (HEJ and the EoS are 6 equa- 
tions, and can determine the 6 parameters of the final 
state, e, n, P, and v. 

Later Csernai [H, 0] pointed out the importance of sat- 
isfying energy, momentum and particle charge conserva- 
tion laws across such hypersurfaces and generalized the 
earlier description of Taub Q to spacelike and timelike 
hypersurfaces (with spacelike and timelike normals re- 
spectively). In this situation the matter both before and 
after the shock was near to thermal equilibrium, and thus 
the conservation laws led to scalar equations connect- 
ing thermodynamical parameters of the two stages of the 
matter: the generalized Rayleigh line and Taub adiabat 

f = [P\(da»d*J/[X] , [P] = [(e + P)X]/(X 1+ X ) . 

(5) 

At much higher energies, at the first stages of the col- 
lision, the matter becomes 'transparent' and the initial 
state is very far from thermal equilibrium. For this stage 
other models were needed to handle the initial develop- 
ment, e.g. refs. Q. The initial non-equilibrium state 
in this situation cannot be characterized by thermody- 
namical parameters or an EoS, so the previous approach, 
with the generalized Rayleigh line and Taub adiabat is not 
applicable. Nevertheless, the intermediate (fluid dynam- 
ical) stage is in equilibrium and has an EoS, while the 
initial state has a well defined energy momentum tensor. 
In this work we will demonstrate that the final invariant 
scalar, thermodynamical parameters can be determined 
in this situation also from the conservation laws. 

Then, Bugaev [l(| E| observed that FO across hyper- 
surfaces with spacelike normals, has problems with neg- 
ative contributions in the Cooper-Frye evaluation [3| of 
particle spectra, thus the FO must yield an anisotropic 
distribution, which he could approximate with a cut- 
Jiittner distribution [Tol . [lT| . This is not surprising as 
in the rest frame of the front (RFF) all post FO particles 
must move "outwards", i.e. p^da^ > is required. This 
condition is not satisfied by any non-interacting thermal 
equilibrium distribution, which extend to infinity in all 
directions even if they are boosted in the RFF. [29| 

Subsequently, another analytic form was proposed by 
Csernai and Tamosiunas, the cancelling- Jiittner distribu- 
tion (l2l | , which replaced the sharp cutoff by a continuous 
cutoff, based on kinetic model results. 

Parallel to this development, the FO process was anal- 
ysed in kinetic, transport approaches |13l4l7| . where the 
FO happened in an outer layer of the spacetime, or in 
principle it could be extended to the whole fluid (al- 
though, at early moments of a collision/explosion, from 
the center of the reaction few particles can escape) . These 
transport studies also indicated that the post FO dis- 
tributions may become anisotropic p^4l6j even for FO 
hypersurfaces with timelike normal [in short: timelike 
surface], if the normal, c&t m , and the velocity four- vector, 
it* 4 , are (very) different. 



These studies led to another FO description, where 
the initial stages of the collision with strongly interact- 
ing matter were described by fluid dynamics, while the 
final, outer spacetime domain (or later times) was de- 
scribed by weakly interacting particle (and string) trans- 
port models, where the final FO was inherently included, 
as each particle was tracked, until its last interaction. It 
is important to mention, that in these approaches, the 
transition from the FD stage to the molecular dynamics 
(MD) or cascade stage happens when the matter crosses a 
spacetime hypersurface, thus the conservations laws 
have to be satisfied and the post FO particle phase space 
distributions [Tol . [l2[ have to be used when the post FO 
distributions become anisotropic. 

In this work for the first time we present a simple co- 
variant solution for the transition problem and conser- 
vation laws for the situations when the matter after the 
front is in thermal equilibrium (i.e. it has isotropic phase 
space distribution) and has an EoS, but the matter before 
the front must not be in an equilibrium state. 

Then we discuss the situation where microscopic mod- 
els are appended to the fluid dynamical model, which are 
in, or close to thermal equilibrium, but the EoS, is not 
necessarily known. 

Subsequently, we present the way to generalize the 
problem to anisotropic matter in final state, which is 
necessary for FO across spacelike surfaces and also for 
timclike surfaces if the flow velocity is large in the rest 
frame of the front (RFF). This problem was solved in 
kinetic approach for the Bugaev cut- Jiittner approach 
[III [l3j and the Csernai- Tamosiunas cancelling- Jiittner 
approach, fl2| by calculating the energy momentum ten- 
sors explicitly from the anisotropic phase space distribu- 
tions, but no general solution is given for post FO matter 
with anisotropic pressure tensor. 



II. NUMERICAL EXTRACTION OF A FREEZE 
OUT HYPERSURFACE 

The transition hypersurface between two stages of a 
dynamical development are most frequently postulated, 
governed by the requirement of simplicity. Thus, such 
a hypersurface is frequently chosen as a fixed coordi- 
nate time in a descartian frame t, or at a fixed proper 
time r from a spacetime point, although in a general 
3+1 dimensional system the choice of such a point is 
not uniquely defined. It is important that the transi- 
tion hypersurface should be continuous, (without holes 
where conserved particles or energy or momentum could 
escape through, without being accounted for). To secure 
that one quantity (e.g. baryon charge) does not escape 
through the holes of a hypersurface is not sufficient, as 
other quantities may (e.g. momentum in case if PdV is 
different on the two sides of a hole). Again, to construct 
such a continuous hypersurface in a general 3+1 dimen- 
sional system is a rather complex task, although, in 1+1 
or 2+1 dimensions it seems to be easy. 
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Both the initial state models and the intermediate 
stage, fluid dynamical models may be such that the cal- 
culation could be continued beyond the point where a 
transition takes place. Then spacetime location of the 
transition to the next stage can or should be decided, 
based on a physical condition or requirement, which may 
be external to the development itself. As a consequence, 
in some cases the determination of transition surface may 
be an iterative process. 

Numerically, the extraction of a Freeze Out (FO) hy- 
persurface is by no means trivial. One of us, BRS, has 
recently provided a proper numerical treatment regard- 
ing the extraction of FO hypersurfaces in two (2D), three 
(3D) and four (4D) dimensions 

For instance, in 2D the history, i.e., the temporal evo- 
lution, of a temperature field of a one-dimensional (ID) 
relativistic fluid can be represented by a gray-level image 
(c/., Fig. 1). In the figure, we use the time t and the 
radius r for the temporal and the spatial dimensions, re- 
spectively. Let bright pixels (i.e., picture elements) refer 
to high temperatures and dark ones to low temperatures 
of the fluid. In this example, a 2D freeze-out hypersur- 
face is an iso-therme. 

In Fig. l.a, we also depict the corresponding co- variant 
normal vectors da^x^ 1 ). In 2D, the length of each nor- 
mal vector is equal to the length of each supporting iso- 
contour vector. Each normal vector has its origin at 
the contra- variant center, x^, of a given contra- variant 
iso-contour vector and points to the exterior of the en- 
closed spacetime region. The latter is also indicated in 
Fig. l.b, where we show that a contra- variant normal 
vector da fl (x fl ) can be obtained by reflection of the co- 
variant normal vector da fM (x fl ) at the time axis (dashed 
line) . 

Finally, in Fig. l.c we show the contra- variant FO con- 
tour vectors with their corresponding contra- variant nor- 
mal vectors da^{x^). Not all of these contra- variant nor- 
mal vectors point to the exterior of the enclosed space- 
time region. 

Note that the sign conventions of the normals of the 
transition hypersurface are important, and must be dis- 
cussed, especially if both timelike and spacelike surfaces 
are studied. In fact, only the timelike contra- variant nor- 
mal vectors point outwards, whereas the spacelike contra- 
variant normal vectors point inwards [22| . 

If we know the FO hypersurface and the local mo- 
mentum distribution after the transition the total, mea- 
surable momentum distribution can be evaluated by the 
Cooper- Frye formula [3J. 



III. EQUATIONS FOR PARAMETERS OF 
FINAL MATTER IN EQUILIBRIUM 

Let us define the contra-variant and co- variant surface 
normal four-vectors as 

da^ = (a t ,a x ,a y ,a z ) = (cr t , <r) 
da^~ (a t , — a x , —a y , — a z ) — (a t ,—a), 

where in general da^da^ = ±D 2 , as the surface ele- 
ment can be either spacelike (-) or timelike (+). We 
can also introduce a unit normal to the surface as: 
da^ = da^ / D so that da^da^ = ±1 . Furthermore 
da» = 7ct (1, Sz) = 7<j(1,s) where for timelike 

surfaces 7 2 = 1/(1— s 2 ) and for spacelike surfaces 
7 2 = l/(s 2 — 1). For the frequently used timelike, one- 
dimensional case da^ = 7,7(1, s) = 7,7(1, 0, 0, dt/dz). 

In the general case the conserved energy-momentum 
current crossing the surface element is 

A" = T» v do v = w u^u" da v — P g^ iV da v (6) 

A^ must be continuous across the freeze-out surface, as 
must the baryon current N^da^, 

j = N^do-ft = n u^da^ = n 7 (a t - (v ■ a)) , (7) 

where j is the invariant scalar baryon charge current. 

We assume that the initial state, "0", and its energy 
momentum tensor and baryon current before the front is 
known. We aim for the characteristics of the final state. 
In total there are six unknowns in the equilibrated final 
state, these are v, e, P and n (here we drop the index 
"1" for the final state for shorter notation), however the 
pressure P, a function of e and n, is given by the EoS, 
P = P{e, n). Knowing n and e, the EoS, and the partic- 
ular form of the corresponding equilibrated distribution 
function, the parameters T, and /1, can also be obtained. 

Thus, we have to solve 5 equations: 



Jo 


= nj(a t - 


-(v- 


<?)) 




(8) 


A ot 


= W7 2 (a t 


-(v 


.«?))- 


Po t 


(9) 


Aq x 


= W7 2 (cr t 


-(v 


• S))v x 


-Po x 


(10) 


Aoy 


= W7 2 (cr t 


-(v 


■ 


-Pa y 


(11) 


A 0z 


= W7 2 (a t 


-(v 


• d))v g 


-Po z 


(12) 



The l.h.s. represents quantities of the initial state of 
matter and the corresponding conserved quantities are 
known. Equations (|9IL0[) can be solved for 7 2 in the cal- 
culational frame: 

2 A ot +Pa t 2 A 0x + Pa x 

w (a t - (v ■ a)) w v x (a t - {v ■ a)) 

Using now Eq. (|10lllll2p one obtains v x , and in a similar 
fashion v y and v z 

_ A 0x + Pa x __ A 0y + Pa y _ A Qz + Pa z 

Vx ~ A Qt + Pa t ' Uy ~ A ot + P<Jt ' Vz ~ A ot + Pa t ' 

(14) 
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FIG. 1: (a) A gray-level image representing the temperature evolution in spacetime of a ID relativistic fluid (see text) which 
is superimposed with contra- variant FO contour vectors (black) and with corresponding co- variant normal vectors da fJ ,(x IJ ') 
(white); (b) co- variant and contra- variant normal vectors, da ll (x ,J ') and d<T M (x M ), (gray) respectively, which originate at the 
contra- variant center x^ (gray) of a single contra- variant FO contour vector (black); (c) as in (a), but with contra- variant 
normal vectors dcr tl (x^) (white) instead of co- variant ones. 



This results for j 2 = 1/(1 — v 2 
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(Apt + P ot) 
{A^ + Pda^) 2 



(15) 



where, (A^ + Pda^) 2 = (A Qt + P a t f - {A 0x + P a x f - 
(Ao y + P <r y ) 2 — (Aq z + P a z ) 2 , is an invariant scalar, 
and 7 transforms as the 0-th component of the 4-vector 
Aq + P G?er M . Notice that eq. ((8]) was not used up to 
this point, thus we can use there results both for the 
baryon-free and baryon-rich case. 

We can have an elegant direct solution for the proper 
energy density, e, and pressure, P, as both of these quan- 
tities are invariant scalars, and we can express these by 
the covariant, 4-vector equation (|6]). From this 4-vector 
equation we can get two invariant scalar equations by (i) 
taking its norm, A^A^^, and (ii) taking its projection to 
the normal direction, A^da^: 

A%A 0li = w 2 (u»da^) 2 + P 2 {da»da^) 
- 2Pw(u' i da tl )u fi g ft "da v 
= w(e - P)(u"da f ,) 2 + P 2 (da ll da tl ) , (16) 

w^da^f - P{d(J^d(J^) . (17) 



A^da^ 



Now expressing w{u^da^) 2 from eq. (I17|) and inserting 
it to eq. (|16p , we obtain our final equation 



(18) 



A$A 0fl = (e - P)A%da ll + e P {da^da^) , 



which can be solved straightforwardly if the EoS, P = 
P(n, e), is known. The other three elements of the equa- 
tion, AqAq^, Agda^, and da^da^, are known from the 
normal to the surface and from energy-momentum cur- 
rent from the pre-transition side. 

Then, eqs. (|13H15[) can be used to determine the fi- 
nal flow velocity. At the end, after all conservation law 



equations are solved, we have to check the non-decreasing 
entropy condition ([3]) to see whether the solution is phys- 
ically possible. If the overall entropy is decreasing after 
transition that would mean that the hypersurface is cho- 
sen incorrectly. One will need to choose more realistic 
condition for the transition and repeat the calculations. 

This result can be used both if the initial state is in 
equilibrium and if it is not. 



A. Final Matter with zero Baryon charge 

In case of an ideal gas of massless particles after the 
front, with an EoS of P — e/3, eq. (TT51) leads to a 
quadratic equation, 



da^da^ e +2 a^da^ < 



3 a^a„ = . 



where a M = Aq/D, is the energy momentum transfer 4- 
vector across a unit hypersurface element. 

If the flow velocity is normal to the FO hypersurface, 
= da^, then for an initial perfect fluid in the Local 
Rest (LR) frame the above covariant equation takes a 
simple form, 



2 eo e — 3 e 







This has two real roots, e = eo (energy density is con- 
served) and e = — 3eo which does not correspond to a 
physical solution, as the energy density should not be 
negative. 



B. Final Matter with Finite Baryon Charge 

If the EoS depends on the conserved baryon charge 
density also, then we must exploit in addition eq. ([7]): 



jo = 3= n^dcrfj) 
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and inserting u^da^ = j/n from here to eq. yields 

j 2 ^ = A%da fl + P{da»da fl ) , 

where w/n 2 = X is the generalized specific volume, well 
known from relativistic shock and detonation theory Q. 
This equation provides another equation for e + P as 

^JL = \ + Piston)] , (19) 

which, together with eq. (IT8|) and the EoS, P = P(e, n), 
provide three equations to be solved for e, P and n. 

This evaluation of the post FO configuration is in 
agreement with the theory of relativistic shocks and det- 
onations |6|, [8| allowing for both spacelike and timelike 
FO hypersurfaces. See also pj. This method of evalu- 
ation observables is frequently used at the end of fluid 
dynamical model calculations (see e.g. (23 - [26| ). 

IV. TRANSITION TO MOLECULAR 
DYNAMICS BEFORE FREEZE OUT 

Recently a frequently practiced method to describe the 
final stages of a reaction is to switch the FD model over to 
a Molecular Dynamics (MD) description at a transition 
hypersurface. This is frequently a fixed time, t, or fixed 
proper time, r hypersurface. The generation of the ini- 
tial state of such an MD model is a task, which depends 
on the constituents of the matter described by the MD 
model. Nevertheless, same principles must be satisfied, 
like the conservation laws, eqs.(Q][2]). 

A. Equilibrium and EoS known before and after 
the transition 

Let us assume, although not required by physical laws, 
that we have thermal equilibrium on both sides of the 
transition and we know explicitly the corresponding final 
momentum distribution of particles. Then, the funda- 
mental equation to construct the post transition micro- 
scopic state, in addition to the conservation laws is the 
Cooper-Frye formula, 

E^=Jfi(x,p)p> t da lt> (20) 

assuming that the local phase space distribution, fi{x,p), 
is known for all initial components of the MD model. If 
fi(x,p) are local equilibrium distributions then (in princi- 
ple) we know the intensive and extensive thermodynam- 
ical parameters and the EoS of the matter when the MD 
model simulation starts. These must not be the same as 
the ones before the transition hypersurface. 

In the usual transition from FD to MD models, where 
the initial state of MD is in equilibrium, the EoS-s are 
known on both sides of the transition surface, and thus, 



both the equations of Ray leigh- line and Taub-adiabat, 
eqs. ([5]), as well as the invariant scalar equations de- 
rived here, eqs. (|16|17|18|19|) can be used to determine 
all parameters of the matter starting the MD simula- 
tion. These then determine the phase space distribu- 
tions, fi(x,p) of all components of the MD simulation. 
Subsequently eq. (|20[) can be used to generate randomly 
the initial constituents of the MD simulation. 

As eq. (|20[) is a covariant equation applicable in any 
frame of reference, the most straightforward is to per- 
form the generation of particles in the calculational frame 
of the MD model. This transition is by now performed 
in many hybrid models combining fluid dynamics with 
microscopic transport models [23]. These models at 
present are the most effective to describe experimental 
data and make the need for a Modified Boltzmann Trans- 
port Equation [la ] less problematic. 

In some cases the first step of the transition, the deter- 
mination of the parameters of the final state from the ex- 
act conservation laws, is dropped with the argument that 
both before and after the transition the matter has the 
same constituents and the same EoS, thus the all exten- 
sive and intensive thermodynamical parameters as well 
as the flow velocity must remain the same. Then, using 
the intensive parameters the final particle distributions in 
the Cooper-Frye formula, eq.(f20|). can be directly evalu- 
ated in a straightforward way. This procedure is correct, 
but only if all features of the two states of the matter and 
their EoS are identical. In some cases the pre transition 
EoS assumes effective hadron masses depending on the 
matter density, while the final EoS is that of a hadron 
ideal gas mixture, but with fixed vacuum masses. This 
leads to a difference in the EoS, thus the above procedure 
is approximate. In such cases, the method can be used, 
but the accurate conservation laws can be enforced by a 
final adjustment step described in the next subsection. 

The situation is similar if the constituents and the EoS 
are almost identical before and after the transition, but 
before the transition a weak or weakening mean field po- 
tential or compression energy is taken into account. 



B. Enforcing Conservation Laws with Approximate 
Generation of the Final State 

In addition to the above mentioned approximate meth- 
ods, even for really identical EoS-s across the transition 
or with generating the final EoS parameter based on con- 
servation laws for the final EoS, inaccuracies may arise 
due to other reasons: during the random generation of 
the initial constituent particles of the MD simulation, the 
exact conservation laws may be violated, due to finite 
number effects. However, the energy and particle num- 
ber conservations are usually enforced during the ran- 
dom generation of particles, even if the above procedure 
of solving the conservation laws beforehand is not fully 
followed. This is usually the consequence of the fact that 
the EoS of the MD model is not necessarily known if the 
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model has complex constituents and laws of motion. 

In any case to remedy this random error and make the 
conservation laws exactly satisfied a final correction step 
is advisable, and it is not always performed. If the energy 
and particle number conservations are enforced then, the 
last variable to balance is the momentum conservation. 
This regulates the flow velocity of the matter after the 
transition initiating the MD simulation. 

The energy momentum tensor and baryon current for 
the generated random set of particle species, "i\ for each 
fluid cell (or group of cells if the multiplicity in a single 
cell is too low) can be calculated from the kinetic defini- 
tion: 

^» = E/ ^rtfpUiM, 
N "(*) = E / ^rffiM- 

which, yield the resulting momentum and flow velocity 
of the matter. This can be used to adjust the flow ve- 
locity to achieve exact conservation of momentum, and 
modify the velocity of generated particles by the required 
Lorentz boost. The other conserved quantities may then 
be affected also, but an iterative procedure to eliminate 
the error completely is not crucial as the error can be 
given quantitatively. 

If the randomly generated state is not following a ther- 
mal equilibrium phase space distribution, fi(x,pi), and 
thus does not have an EoS, the above described scalar 
equations cannot be used to generate the initial config- 
uration of the MD model. Nevertheless, the second step 
to check the conservation laws with the kinetic definition, 
and then correct the parameters of the generated parti- 
cles can be done. For a required level of accuracy in this 
case an iterative procedure may be necessary. 

Another, easier way to remedy this problem is to 
choose the transition hypersurface earlier so that the sub- 
sequent matter is still in thermal equilibrium. This can 
always be done if the requirement of entropy increase is 
satisfied. 



V. FINAL STATE OUT OF THERMAL 
EQUILIBRIUM 

We have mentioned that the assumption for having 
thermal equilibrium in the final state is neither excluded 
nor required from transport theoretical considerations. 
However, thermal equilibrium distribution is not possible 
if we have to describe FO across a spacelike hypersurface 
(see the discussion in section HJ) 

In the MD model description the final post FO momen- 
tum distributions develope a local anisotropy if the FO 
has locally a preferred direction. Unless the unit normal 
of the FO hypersurface is equal to the local flow velocity 
of the pre FO matter, there is always a selected spatial 



direction which is the dominant direction of FO. This sit- 
uation is discussed in several theoretical works, and some 
general features can be extracted from these studies. 



A. Approximate kinetic models for Freeze Out 

In explicit transport models this situation is handled 
[irj| - [l3| : starting from an equilibrium Jiittner distri- 
bution and considering a momentum dependent escape 
probability in the collision term, - which reflected the di- 
rection of the FO front and the distance from the front, - 
an anisotropic distribution was obtained (i.e. a distribu- 
tion, which was anisotropic even in its own LR frame). 

This anisotropic distribution could be approximated 
with analytic distribution functions The starting 

point is the un-cut, isotropic, Jiittner distribution in the 
rest frame of the gas (RFG), which is centered around 
the 4- velocity vector, u RFG . This distribution is then 
cut or cut and smoothed. The resulting distribution has 
a different new flow velocity, u LR , which is non-zero in 
RFG, and is pointing in space in the direction of the 
normal of the FO hypersurface, a, labeled by ||. This 
u^ LR defines the Local Rest (LR) frame of the post FO 
matter. 

The spatial direction of a is not affected by the Lorentz 
transformation from RFF to RFG and then to LR, as a 
is the direction of the Lorentz transformation from RFG 
to LR [^l . In the general case the boost in the vrfg,± 
direction leads to a change of the distribution function in 
the p± direction, but does not affect the distribution in 
the direction, or the procedure of cutting or cancelling 
the distribution in the || direction. (The illustration in 
Fig. 2a shows the spatial momentum distribution where 
the boost in the orthogonal direction vrfg,±- is already 
performed.) 

In the final LR frame, the matter is characterized by 
a rather complex energy momentum tensor, inheriting 
some parameters from the original uncut distribution in 
RFG, like the temperature and chemical potential, but 
as the resulting distribution is not a thermal equilibrium 
distribution, these parameters are not playing any ther- 
modynamical role. One has to determine all parameters 
numerically from conservation laws (III2I) . as done in refs. 
01. 

Interestingly, a simplified numerical kinetic FO model 
13] led to a FO distribution satisfying the condition 
p^da^ > for spacelike FO with a smooth distribution 
function, which is anisotropic (also in its own LR frame) 
and has a symmetry axis pointing in the dominant FO 
direction. This distribution was then approximated with 
an analytic, "cancelling- Jiittner" distribution which 
can also be used to solve the FO problem. 

After FO, the symmetry properties of the energy mo- 
mentum tensor are the same for the cut- Jiittner and 
cancelling- Jiittner cases The FO leads to an 

anisotropic momentum distribution and therefore to an 
anisotropic pressure tensor. The energy momentum ten- 
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sor is not diagonal in the RFG frame, there is a non van- 
ishing transport term, T° l [lj [l2[ , in the 2-dimensional 
plane spanned by the 4- vectors, u RFG and da^. One 
can, however, diagonalize the energy momentum tensor 
by making a Lorentz boost into the LR frame using Lan- 
dau's definition for the 4- velocity, u^ R . In this frame then 
the energy momentum tensor becomes diagonal, but the 
pressure terms are not identical, due to the anisotropy of 
the distribution: 



T» v = diag(e,P,|,Px,P± 



(21) 



LR 



Here the energy density, e, of course must not be the 
same as in the case of an isotropic, thermal equilibrium 
post FO momentum distribution. This can be seen from 
the kinetic definition of the energy momentum tensor as 
shown in refs. [TlL [T^| . 

We need the complete post FO momentum distribu- 
tion and the corresponding energy momentum tensor to 
determine final observables. This depends on the trans- 
port processes at FO, and cannot be given in general; 
however, due to the symmetries of the collision integral, 
the symmetries of the energy momentum tensor are the 
same irrespectively of the ansatz used (e.g. cut-Juttner, 
cancelling- Jiittner or some other distribution). 

In kinetic transport approaches the microscopic escape 
probability 15] is peaking in the direction of d&^, which 
yields a distribution peaking in this direction, i.e. yield- 
ing the same symmetry properties as the previously men- 
tioned analytic ansatzes. The energy momentum tensor 
in general takes the form 



LR LR 



LR, 



(Pi - P ± )F»F V (22) 



where A^, is the orthogonal projector to u LRl and P^ 
is the unit 4- vector projection of da^ in the direction or- 
thogonal to u1 R , i.e. P M = CA^ R da u , where C ensures 
normalization to -1. In the Landau LR frame this returns 
expression ([2~Tjl . The 4- velocity, u^ R , and the other pa- 
rameters of the post FO state of matter, should be deter- 
mined from the conservation laws (|1I2|) . The schematic 
diagram of the asymmetric distributions and the different 
reference frames can be seen in Fig. [5] 

The FO problem was solved for these configurations 
and ansatzes, by satisfying the conservation laws explic- 
itly for the full energy momentum tensor. We do not have 
a general EoS(s) that would characterize the connection 
among e, Py and P±, furthermore the relation connecting 
these quantities depends on the 4-vectors da^ and u^ R . 
In addition this connection depends on the details or as- 
sumptions of the transport model. The simple models 
[Til fl2j provide examples for such a dependence. If da^ 
is known, then for baryon free matter we can determine 
four unknowns: u^ R and, an additional parameter of the 
post FO distribution, from eq. (Due to normal- 

ization only 3 components of u LR are unknowns.) For 
baryon-rich matter we can determine one more unknown 
parameter, since we have one additional equation, the 
conservation of baryon charge from eq. (H} . 



B. Exploiting general symmetries of anisotropic 
final states 



The first step of solution can be done similarly to the 
isotropic case. Then in eq. © the enthalpy will change 
as id -> el Pi = ioj_j and P — > Pj_, plus an additive 
term will appear, (P|| — P±) F^F 1 ' da v . Furthermore, 
eqs. ([9l ll2p remain of the same form, with w±, and Pj_, 
plus the additive term (Pi — PjJ P^P" da v will appear 
in the r.h.s. of cqs.(|9l [T"2"]) . This additive term will also 
appear in the expression of v x after eq. (I13j) and in the 
denominator of eq. (fl~5|) also. 

The additional term, (P|| - P ± ) F^F" da v , in eq. © 

is orthogonal to u M (by definition of P M ), so when we cal- 
culate the scalar product (fTBjl their cross term vanishes, 
so 



4Mo„ = 



A^da^ 



w x (e - P ± ){u»da^) 2 + (P ± ) 2 (da»da^ 
(P,l -P±)(P|| + P X )(F» da v ) 2 , (23) 



(P|,-Px)(P i/ da v f 



(24) 



Now one can express w±(u ,J 'da f j l ) 2 from eq. ([24| and in- 
serting it to eq. (f23|) . we obtain that 

A%A „ = (e - P±)Agdtrp + e P ± {da^da^) 

-(e + P ] \)(P\\-P ± )(F"' da,) 2 . (25) 

where this equation is not a scalar equation as it dependcs 
on P M = CA^dtTj,, where the projector is dependent on 
w' 1 . These equations are similar to the ones obtained for 
the isotropic case, however, to solve this last equation we 
need a more complex relation among e, Pi, P±. As these 
arise from the collision integral in the BTE approach the 
needed relation may depend on Uq and da^. On the 
other hand, the escape probability may be simple, or 
may be approximated in a way, which yields an ansatz 
for this relation with adjustable parameters, and then the 
problem is solvable. This was the case in refs. 0,(3. 

The recent covariant formulation of the kinetic freeze- 
out description [15] indicates that the relation among the 
different parameters of the anisotropic energy momen- 
tum tensor, should be possible to express in terms of 
invariant scalars, which may facilitate the solution of the 
anisotropic FO problem. 

When the adjustable parameters of the post FO matter 
are determined in this way from the conservation laws, 
we still need the underlying anisotropic momentum dis- 
tribution of the emitted particles in order to evaluate the 
final particle spectra using the Cooper- Frye formula with 
this anisotropic distribution function. Once again, when 
all conservation law equations are solved we have to check 
the non-decreasing entropy condition to sec whether the 
selected FO hypersurface is realistic. 

In case of an anisotropic final state, due to the in- 
creased number of parameters and their more involved 
relations, the covariant treatment of the problem may 
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FIG. 2: (Color online) (a) Illustrative contour plot of momentum distributions of particles after freeze-out (FO) to an anisotropic 
final state. Here the py-axis points into the direction of the space-like normal of the FO hypersurface, while pj_ is the orthogonal 
direction to p«. Line "1" indicates the un-cut Jiittner distribution (dotted-line) , "2" indicates the cut-Jiittner distribution 
(dashed-line), "3" indicates the canceling- Jiittner distribution (dashed-dotted line), "i?" is the center of the spherical un-cut 
Jiittner distribution moving with velocity, u RFG , and "L" is the center of the cut and canceling- Jiittner distributions, which 
move to the right along the parallel momentum direction with velocity, u^ R , corresponding to their LR frame. If u RFG and u^ LR 
have a non-vanishing _L component the resulting distribution can be obtained by an additional Lorentz boost in that direction, 
(b) Space-time reference frames for anisotropic final state distributions. The time axis, t, and the spatial axis, x«, represent 
the "Rest Frame of Front" (RFF), which is the defining front for the cut-Jiittner and canceling- Jiittner distributions. "L" 
indicates the "Local Rest Frame" (LR) of the post FO cut or canceling- Jiittner distribution, "i?" indicates the "Rest Frame of 
Gas" (RFG) corresponding to the isotropic, un-cut Jiittner distribution. Thus the velocities u RFF , u RFG and u^ R point in the 
directions of the time axes of the corresponding reference frames, indicated by t, t , t respectively. d<r M is the normal vector 
of the front, it has the same direction as the spatial axis, x«. The spatial axes of the frames RFG (R) and LR (L) are shown 
with the same line style (color) as the corresponding distributions. The LR frames for case 2 and 3 are the same. 



not provide a simplification, compared to the direct so- 
lution of conservation laws for each component of the 
energy momentum tensor (e.g. (TTI. Il2j). 

C. Anisotropic initial and final states 

Recent viscous fluid dynamical calculations evaluate 
the anisotropy of the momentum distribution is in the 
pre FO viscous flow (see e.g. [28[ .) This anisotropy is 
governed by the spacetime direction of the viscous trans- 
port. The pre and post FO matter may still be different, 
e.g. the pre FO state may be viscous QGP with cur- 
rent quarks and perturbative vacuum, while post FO we 
may have a hadron gas or constituent quark gas. The 
final state will also be anisotropic, not only because of 
the initial anisotropy but also due to freeze-out. The two 
physical processes leading to anisotropy are independent, 
so their dominant directions are in general different. In 
this case the general symmetries are uncorrelated and 



cannot be exploited to simplify the description of the 
transition. Due to the change of the matter properties, 
the conservation laws, eqs. (JTJ3JI, are needed to deter- 
mine the parameters of the post FO matter before the 
Cooper-Frye formula with non-equilibrium post FO dis- 
tribution is applied to evaluate observables. 

VI. SUMMARY 

In this work a new simple covariant treatment is pre- 
sented for solving the conservation laws across a transi- 
tion hypersurface. This leads to a significant simplifica- 
tion of the calculation if both the initial and final states 
are in thermal equilibrium. The same method can also 
be used for the more complicated anisotropic final state, 
however, this method is only advantageous if the more 
involved relations among the parameters of the post FO 
distribution and the distribution itself is given in covari- 
ant form, preferably through invariant scalars. 
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